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We introduce cluster dynamical models of conflicts in which only the largest cluster can be involved 
in an action. This mimics the situations in which an attack is planned by a central body, and the 
largest attack force is used. We study the model in its annealed random graph version, on a fixed 
network, and on a network evolving through the actions. The sizes of actions are distributed with 
a power-law tail, however, the exponent is non-universal and depends on the frequency of actions 
and sparseness of the available connections between units. Allowing the network reconstruction 
over time in a self-organized manner, e.g., by adding the links based on previous liaisons between 
units, we find that the power-law exponent depends on the evolution time of the network. Its lower 
limit is given by the universal value 5/2, derived analytically for the case of random fragmentation 
processes. In the temporal patterns behind the size of actions we find long-range correlations in the 
time series of number of clusters and non-trivial distribution of time that a unit waits between two 
actions. In the case of an evolving network the distribution develops a power-law tail, indicating 
that through the repeated actions, the system develops internal structure which is not just more 
effective in terms of the size of events, but also has a full hierarchy of units. 

PACS numbers: 89.75.-k; 05.65. +b; 89.20.-a 



I. INTRODUCTION 

The study of coagulation-fragmentation models in 
physics have a long history, and have been used to ex- 
plain a variety of physical phenomena, including the for- 
mation of aerosols, colloidal aggregates, polymers and 
celestial bodies More recently the basic theory, orig- 
inally developed by Smoluchowski [2| and Becker-Doring 

has been adapted to model group behaviour in human- 
dynamics. The first microscopic model of herding in fi- 
nancial markets consisted of randomly connected agents 
[3| • In this model agents are connected with probability a 
and disconnected with probability 1 — a. Agents that are 
in the same group share information and make the same 
decisions in the market. The parameter a was tuned to 
the percolation threshold to obtain a power-law distribu- 
tion of group sizes and hence a power-law distribution of 
returns. 

Inspired by this, in 2000 Eguiluz and Zimmermann 
[| (EZ) introduced a coagulation-fragmentation model 
of herding in financial markets. The EZ model is a ki- 
netic version of the model in [J] , in which the group sizes 
emerge more naturally in the limit of large time. At each 
time step with probability p an edge is introduced be- 
tween two randomly selected agents or with probability 

1 — p the group of a randomly selected agent is frag- 
mented. In Dhulst and Rodgers solved this model 
exactly and showed that in the limit t — >■ oo the system 
evolves to a stationary cluster size distribution with the 
number of clusters of size s, n^, is given by 



p'*-i(2s-2)! 



N 



(1) 



where N is the number of agents in the system. In the 



limit p — 1 this result can be expanded for large s to 
give a power-law with Us ~ s~'^ and r = 5/2. 

This basic model has been used Q to explain the ap- 
parent ubiquity of the exponent r = 5/2 in the data from 
modern insurgent warfare. Many authors have studied 
wars and confiicts empirically, and many have concluded 
that the distribution of daily casualties, or of casualties 
per attack, has a power-law distribution [1, A wide 
range of exponents have been reported. In [10[ an ex- 
ponent of r = 1.8 is reported for the intensity of 119 
old wars between 1816 and 1980. Casualty numbers in 
global terrorist attacks, since 1968, have r = L7 for G7 
countries and r = 2.5 for non-G7 countries [il|. More 
recently, in [7i] the daily data from killings and injuries in 
Columbia and civilian casualties in Iraq were examined. 
The total data sets gave good power-law distributions 
with r = 3 for Columbia and r = 2 for Iraq. By reducing 
the time window within these data sets, and then sliding 
the time window forward, the authors of fz^ were able 
to calculate a series for r as a function of time. This 
revealed that as time window is shifted towards recent 
time the exponent r in both the Columbian and Iraqi 
datasets was tending to the value r = 5/2. This caused 
the authors of Q to develop a new model of modern in- 
surgent warfare based on a dynamical model of herding 
introduced in d, @ . 

In financial markets, the return is equal to the differ- 
ence in the number of buyers and sellers at a particular 
time, or alternatively to the supply and demand balance 
at a particular time [ll] . Hence in the models of finan- 
cial markets the probability of having a return of size s 
is proportional to sus. In models of wars and conflicts 
the units can be thought of as troops, weapons or equip- 
ment and it is an assumption in these models ffl that 
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the number of casualties inflicted by a force of strength s 
is proportional to s. In these models the parameter p is 
an external fixed parameter, however, in the real dynam- 
ics this probability might emerge from another stochastic 
process, inherently related to the main dynamics. One 
would expect that the stochastic processes would be quite 
different in the market and conflict dynamical systems. 

The empirical studies of war and terrorist actions 
demonstrate that the power-law behaviour is non- 

universal, with the exponent r depending on the type 
of conflict (insurgence, guerrilla, terrorist, etc..) and 
on the geographical location and period of time consid- 
ered. In contrast the original model with the coagulation- 
fragmentation of random clusters always leads to a uni- 
versal exponent r = 5/2. Thus the EZ model is not 
able to capture the underlying mechanisms which lead 
to the variety of observed power-laws. Moreover, the full 
fragmentation of clusters, which is suitable for herding 
in market dynamics, might not be appropriate for the 
conflicts. 

In this paper we study new models of conflicts with 
the cluster dynamics on networks when only the largest 
cluster can be involved in an action. This mimics the sit- 
uations in which an attack is planned by a central body, 
and the largest attack force is used, and where the units 
interact within a social or technological network. Follow- 
ing the action, the network is either reset back to the 
original structure (fixed network model), or evolves by 
adding new links between units involved in the action 
(emergent network model). The evolving network model 
provides a mechanism for understanding the power-law 
exponent dependence on the evolution time, similar to 
the empirical data. Furthermore, in our simulations we 
analyse the temporal behavior of each unit, which is be- 
hind the emergent power-laws. This allows us to unravel 
the changing nature of the group actions and monitor 
emergent structure of their connections. For comparison, 
we also study the cluster dynamics with maximal cluster 
fragmentation on an annealed random graph structure, 
which is closer to the original model and the situation in 
real market dynamics. 

In Section [ll] we introduce the microscopic rules of 
the model and present the simulations on the annealed 
random graph, the situation which parallels the origi- 
nal coagulation-fragmentation model of market crashes. 
We also define and compute the temporal quantities that 
can be monitored in real data. In section IIIII we present 
simulation results corresponding to the cluster dynamics 
on a network with fixed and evolving network topology. 
Finally, in section IIVI we give a brief summary of the re- 
sults and a discussion, in particular on the mechanisms 
revealed within our model, which might drive the evolu- 
tion of modern conflicts. 

II. CLUSTER-AGGREGATION AND LARGEST 
CLUSTER FRAGMENTATION WITH AN 
ANNEALED RANDOM GRAPH TOPOLOGY 

In this section we consider a minimal extension of the 
original aggregation-fragmentation model @| , in which 



instead of choosing a random group to fragment we al- 
ways choose the largest group. At each time step with 
probability p a pair of units is selected and joined to- 
gether, whereas with probability 1 — p the largest cluster 
is found and fragmented. As in the original model of 
market dynamics 0, @, in the system of N units any 
randomly selected pair can join each other in the aggre- 
gation event, whereas all the links between pairs within 
the cluster are removed when the cluster is destroyed. In 
terms of topology (see later) this situation is represented 
by a random graph with an annealed link structure. The 
effects of a fixed and evolving network underlying the 
aggregation processes will be studied in section IIIII 



A. Cluster size statistics and temporal structure 

Starting with TV = 1000 units, we apply the rules 
of random aggregation and fragmentation events, as de- 
scribed above. If we always select the largest cluster for 
fragmentation, instead of a randomly selected cluster, it 
affects the system dynamics in a different ways, depend- 
ing on potential size that the cluster can grow, which, 
on the other hand is determined by the fragmentation 
probability 1 — p. For a sizable probability of fragmenta- 
tion, the system has no large clusters, which causes the 
probability distribution of cluster sizes to decay rapidly. 
The distribution of all sizes of clusters and the size of the 
fragmented cluster is shown in Fig. [TJi for several values 
of the probability 1 - p = 10%, 5%, and 2%. It exhibits 
a power-law tail for sizes above a threshold value s > sq 
with the scaling exponent t defined by the expression 

P{s)^s-\s>SQ. (2) 



In each case, the size of the largest (fragmented) cluster 
starts playing a dominant role over a threshold sq (max- 
imum of the fragmented size distribution), which moves 
upward with reduced fragmentation probability. The 
slope T of the cluster-size distribution over the threshold 
is fully determined by the distribution of the fragmented 
cluster size and it is continuously varying with the proba- 
bility p. With the reduced fragmentation probability (in- 
creasing aggregation probability p) the slope of the size 
distribution reduces and converges towards the universal 
curve obtained in the random cluster fragmentation (top 
line in Fig.[TJi). Intuitively, this is understandable since, 
for vanishingiy small but finite probability of crush, a 
very large clusters may occur in the system increasing 
the probability to be selected by random picking unit. 
However, the nature of correlations that the system de- 
velops over time reveals that these are different processes. 
The time fluctuations in the number of clusters and their 
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FIG. 1: For the annealed random-graph topology: Time series (c) and their power spectra (a) of fluctuating number of 
clusters for aggregation probability p = 0.9 and two crush mechanisms — maximal cluster crush and random cluster crush, 
(d) Distribution P(s) of size s of all clusters, empty symbols, and size of crushed maximal clusters, crosses, for three values 
of parameter p = 0.9,0.95, and 0.98. Dashed lines indicate slopes 6, 4, and 3, of the corresponding curves. Also shown is 
the distribution of size in the case of the destruction of random clusters (o), with full line indicating the slope 2.5. (b) The 
distribution P{At) of waiting time At between fragmentation events, shown with full symbols, and the distribution of waiting 
times of individual units involved in the events, empty symbols, for same values of p and the fragmentation of the maximal 
cluster. Large data are logarithmically binned. 



power-spectra are shown in Fig.[T^,c for fixed aggregation 
probability. Apart from the visual difference, in partic- 
ular in the sizes of destruction events and the average 
number of clusters present (697 in random, 628 in max- 
imal cluster destruction events), the power spectrum of 
these time series shows long-range correlations with dif- 
ferent exponents. Specifically, the power spectrum S^i^) 
decays at high frequencies ly as 

sH ^ (3) 

with the scaling exponent (f> depending on the fragmen- 
tation mechanism: We find the exponent close to (j) 1.6 
for the maximal cluster fragmentation, compared to the 
short-range correlations in the case of the random cluster 
fragmentations 1.9 and same aggregation probability 
p = 0.9. 

Further insight into the mechanisms of the cluster dy- 
namics can be gained from the temporal statistics, in 
particular, the distribution of the time intervals between 
events. This distribution can be looked at the level of the 
whole system, as well as the level of each unit. The re- 
sults of the simulations are shown in Fig. \Vp. The distri- 
bution between the fragmentation events is exponential, 
refiecting the fact that in the model the fragmentation 



occurs randomly with the external parameter 1 —p. How- 
ever, from the point of view of the individual units, the 
distribution of the time intervals between two successive 
fragmentation events in which a unit has been involved 
is rather non-trivial. The distribution seems to be pri- 
marily determined by the topology (a random annealed 
graph in this case) where the linking between the units 
occur, rather than with the probability of fragmentation 
(see also discussion in the next section). 



III. CLUSTER DYNAMICS ON NETWORKS 
WITH FIXED AND EMERGENT TOPOLOGY 

In this section we test the effect of placing the aggrega- 
tion process on a network. In particular, the units are as- 
sociated with network nodes, while the network structure 
permits the aggregation only along the links between the 
units. Thus, in contrast to the annealed random graph 
topology discussed above, on the network a cluster can 
grow by aggregating with locally available units. Apart 
from the network structure constraints in the aggregation 
process, the action event occurs with the probability 1 —p 
always involving the largest cluster as above. After each 
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action event, the units who were involved in the action 
can build new connections, specifically, they form a clique 
of the size of the cluster, which then plays a role in future 
aggregation. Thus, the dynamical reconstruction of the 
network structure is tightly linked to the action events. 
We consider the temporal and statistical features of the 
process and the structure of the emergent network. Here 
we present the simulation results starting from a tree net- 
work which contains several subtrees with local hubs (the 
algorithms to grow such networks are described in [l3|). 
An example of the network structure emerging through 
the cluster is shown in Fig. [5] along with the network 
representing the situation of the annealed random graph 
topology, studied above in secHIl As shown below, these 




FIG. 2: Emergent networks in two types of cluster dynamics: 
(top) annealed random graph topology, and (bottom) network 
built via action events on the tree-of-trees structure. 

rules with the aggregation on network and action events 
with the reconstructed links appear to be a suitable basis 
for describing the evolution-dependent power-law expo- 
nent (which was observed for example in the analysis of 
the empirical data of war events in Columbia Q). For 
comparison, we also consider the case when the network 



structure remains fixed throughout the entire process. In 
tis case after the action involving the largest cluster, the 
links between involved units are always reset to the orig- 
inal network structure. 




FIG. 3: Time series (bottom) and power spectrum (top) of 
the number of clusters in the case of maximal-cluster action 
mechanism on the evolving network for two values of aggre- 
gation probability probability p. 

The network reconstruction with building new connec- 
tions between units involved in the action is affecting the 
dynamics through increased temporal correlations. In 
Fig. [3] the example of two time series of number of clusters 
are shown together with their power spectra, correspond- 
ing to two values of the parameter p. We find the power- 
spectrum correlations at hight frequencies with the expo- 
nent (f) = 1.75 ± 0.026 for p = 0.9, while (/) = 1.36 ± 0.035 
for p = 0.95. Further increase of the correlations towards 
the 1 // noise behavior is expected in the asymptotic limit 
of small action probability. In both cases only weak cor- 
relations occur at small frequencies (long times). 

The distribution of cluster sizes, shown in Fig. 21 also 
varies with p with the tail slope dominated by the distri- 
bution of the maximal cluster involved in the action. The 
exponents in the case of the emergent network built on 
the tree-of-trees structure are found to be smaller com- 
pared to the maximal cluster sizes on the annealed ran- 
dom graph. Specifically, the exponents are 6.3, 3.8, and 
3.1 for p =0.9, 0.95 and 0.98, respectively. The results 
are for the evolution time T = 10^ steps. Simulations 
for different evolution times lead to different scaling ex- 
ponents, as shown in the inset to Fig. |4l Within our 
model the temporal dependence of the exponent r is re- 
lated to the evolving network structure. Specifically, the 
network is reconstructed by gaining new links after each 
action event, as explained above, thus the average con- 
nectivity of units increases over time. In the average, the 
number of the action events for the evolution time T is 
given by (1 ~ p)T. The increased connectivity enables 
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larger clusters to form, however, this occurs at the ex- 
pence of smaller clusters. Consequently, the distribution 
is steeper for longer evolution time. The decrease of the 
exponent with decreasing evolution time can be also un- 
derstood by considering the aggregation probability as 
follows: In order to have the same network structure sta- 
tistically when the evolution time T is doubled, T — )■ 2T, 
one needs to reduce the number of the reconstruction 
events per time step, i.e., 1 — p — > (1 — p)/2. This means 
that the aggregation probability is effectively larger and 
thus the exponent closer to 5/2. 
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FIG. 4: Size distribution of all clusters on the evolving net- 
work for p =0.9,0.95, 0.98, and evolution time T = 10^ steps. 
Full lines indicate slopes r =6.3,3.8, 3.1, while dashed line has 
the slope r = 2.5. Inset: Scaling exponent r plotted against 
inverse evolution time T for fixed p = 0.95. 

Further comparative analysis of the network effects on 
the cluster dynamics is given in Fig. [5] Here we also con- 
sider the case where the underlying tree network remains 
fixed in the actions, i.e, following the action events only 
the original links are preserved. It is remarkable that, 
when the largest cluster starts dominating the dynam- 
ics (above the threshold value Sq ~ 6 in this case), the 
tail of the distribution on the tree network coincides with 
the tail obtained in the annealed random graph topology. 
This supports our conclusion that the local structure of 
the network plays a pivotal role in this cluster dynamics. 
Specifically, the tree graph structure and the annealed 
random graph, which is locally tree-like, have the same 
effects statistically on the cluster aggregation and conse- 
quently on the size of the largest cluster. 

Apart from the cluster size distribution, wc find that 
the network structure also plays a role in how often the 
units participate in the action. The distribution of wait- 
ing times between successive actions of a unit, averaged 
over all units on the network, is shown in Fig. [SJtop) 
for the three network topologies. Here we see differences 
appear between annealed random graph and fixed tree 
network. Apart from small times, the data can be fitted 
by the g— exponential function (see [l3| for the origin of 
this distribution in long-range interaction systems): 

P{M) = B(At)"[l + (1 - g)At/ro]i/(i-«) , (4) 



FIG. 5: (top)Distribution of time intervals between consec- 
utive action events of a particular unit, averaged over all 
1000 units, for three network topologies: annealed random 
graph, fixed network of tree-with-subtrees, and evolving net- 
work built on the tree-with-subtrees. (bottom) Distribution 
of sizes of all clusters for the same networks. Fixed aggrega- 
tion probability p = 0.9 and simulation time T = 10'' steps. 

with a prefactor and different parameters. Specifically, 
q =1.04 (exponential limit!) for the random graph topol- 
ogy, and q =1.12 for the fixed tree network. Whereas, 
a prominent power-law tail with the with the slope 
T =2.66, compatible with q =1.38, occurs in the case 
of the evolution of network. This implies that the units 
which are deployed less often than the threshold time 
(To ~ 10'^ steps), appear to be self-organized in a hierar- 
chical manner. 



IV. CONCLUSIONS 

We have studied the cluster aggregation- 
fragmentation processes on networks and assuming 
that the largest available cluster is always involved in 
the fragmentation or action events, which also may 
affect the network structure. This is motivated by the 
dynamics of conflicts, where planned war actions deploy 
the largest forces (and may result in correspondingly 
large damage). The underlying network with its fixed 
or dynamically evolving structure represents constraints 
for the aggregation processes, thus affecting the largest 
cluster sizes, which is relevant for war conflicts. We also 
studied the situation where such constraints do not exist 
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(e.g., in the herding events in stock market dynamics, 
where the information is globally available), resulting 
in an annealed random graph structure. The process is 
controlled by the aggregation probability p, which may 
result from another coupled stochastic process, but here 
it is taken as an external parameter. 

We find a number of new features, in both the geo- 
metric and temporal patterns of the process, which can 
be related to the network topology. In particular, the 
system is characterized by: 

• non-universal power-laws in the distributions of the 
sizes of events; the scaling exponents appear to vary 
with the fragmentation probability and the evolu- 
tion time (in the case of the evolving network); 

• long-range temporal correlations in time series rep- 
resenting fluctuations of the number of clusters; 
the correlations are intensified on networks whose 
structure evolves through the action events; 

• evolution of the internal organization, i.e., connec- 
tivity and hierarchy between units, as a result of 
the action events. 

The distribution of cluster sizes in the mechanism when 
the largest cluster is always selected for the action, is 
conditioned by the local connectivity of the network. Ag- 
gregation process on trees yields large number of small 
clusters due to sparse connections on the tree. Lack of 
large clusters results in the large scaling exponent of the 
size distribution. Within network terminology, the clus- 
ter aggregation-fragmentation processes in "free space" 
leads to a sparse random graph structure with annealed 
links. Its local structure is also tree-like, leading to sta- 
tistically similar cluster size distribution as on the fixed 
tree graphs. On the other hand, on the evolving network 
which gains new links (cliques) after the action events, 
the node connectivity is progressively increasing. This 
allows larger clusters to form, and consequently, smaller 
scaling exponent of the size distribution, compared to 
the sparse networks for the same fragmentation proba- 
bility. In both cases the exponent decreases when the 
fragmentation probability 1 — p is reduced. In the limit 



p — 1 with a vanishingly small but finite fragmentation 
probability, the distribution is expected to approach the 
classical slope 5/2 of random aggregation-fragmentation 
processes, however, the cutoff might be different. This 
limit needs more theoretical consideration, which is left 
out of the present work. The relation between the net- 
work connectivity and the role of the largest cluster in 
the dynamics is also the basis for the observed variation 
of the scaling exponent of all clusters r with the evolu- 
tion time of the network in our model. For a fixed frag- 
mentation probability (number of reconstruction events 
per time step), the larger evolution time leads to larger 
size of the fragmented cluster. Hence the distribution of 
all clusters is steeper, compared to small evolution time, 
where the difference between cluster sizes is reduced and 
thus closer to the processes with random cluster events. 

Our analysis of temporal features of the system reveals 
the nature of the process underlying the observed clus- 
ter size distribution. Namely, systematic network recon- 
struction following the actions over time leads to a com- 
plex organization between units which causes them to 
play different roles (frequencies) in the actions. Despite 
the temporal correlations, which are caused by the in- 
volvement of the largest cluster in the events, such inter- 
nal structure does not occur in the original aggregation- 
fragmentation model in free space and locally sparse net- 
works. Thus the network reconstruction via the clus- 
ter dynamics proposed in our work provides the mech- 
anism for building an effective structure of units, which 
has large capacity (in terms of the size of events) and a 
self-organized internal structure. 
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